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Abstract 

We consider the diffractive Drell-Yan process in proton-(anti)proton collisions at high energies in 
the color dipole approach. The calculations are performed at forward rapidities of the leptonic pair. 
Effect of eikonalization of the universal "bare" dipole-target elastic amplitude in the saturation 
regime takes into account the principal part of the gap survival probability. We present predictions 
for the total and differential cross sections of the single diffractive lepton pair production at RHIC 
and LHC energies. We analyze implications of the QCD factorisation breakdown in the diffractive 
Drell-Yan process, which is caused by a specific interplay of the soft and hard interactions, and 
resulting in rather unusual properties of the corresponding observables. 
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I. INTRODUCTION 



The exclusive diffractive production of particles in hadron-hadron scattering at high ener- 
gies is one of the basic tools for both experimental and theoretical studies of the small- x and 
nonperturbative QCD physics. The understanding of the mechanisms of inelastic diffraction 
came with the pioneering works of Glauber [lj], Feinberg and Pomeranchuk Good and 
Walker [3] . If the incoming plane wave contains components interacting differently with the 
target, the outgoing wave will have a different composition, i.e. besides elastic scattering a 
new diffractive state will be created (for a detailed review on QCD diffraction, see Ref. |4|). 
Among the most important examples, the leading twist diffractive Drell-Yan (DDY) process 
is of a special interest since it gives rise to a clean experimental signature for the QCD 
factorisation breaking effects where soft and hard interactions interplay with each other jBf, 
thus, providing an access to the soft QCD physics 0]. 

Typically, the single-diffractive Drell-Yan reaction in pp collisions is characterized by a 
relatively small momentum transfer between the colliding protons, such that one of them, 
e.g. pi, radiates a hard virtual photon k 2 = M 2 ^> m 2 and hadronizes into a hadronic 
system X both moving in forward direction and separated by a large rapidity gap from the 
second proton P2, which remains intact, i.e. 

Pi+p 2 ^l*(l + l~)+X + (gap)+p 2 (1.1) 

Both the di-lepton and X, the debris of pi, stay in the forward fragmentation region. In this 
case, the virtual photon is predominantly emitted by the valence quarks of the proton p\. 
Below we will refer to this as the diffractive Drell-Yan process at forward rapidities. Notice 
that this is different from double diffractive Drell-Yan process, where the di-lepton l + l~ is 
produced at central rapidities, while both protons survive the collision (see e.g. Ref. 
Then, the 7* can be emitted by a sea quark or antiquark. We postpone this case for future 
studies, and concentrate here on the single diffractive Drell-Yan process. 

In some of previous studies Refs. [7|, |8j of the single diffractive Drell-Yan reaction the 
analysis was made within the phenomenological Pomeron-Pomeron and 7-Pomeron fusion 
mechanisms using the Ingelman-Shlein approach [9[ based on QCD and Regge factorization. 
This led to specific features of the differential cross sections similar to those in diffractive 
DIS process, e.g., a slow increase of the diffractive-to- inclusive DY cross sections ratio with 
c.m.s. energy y/s, its practical independence on the hard scale, the invariant mass of the 
lepton pair squared, M 2 jzj. 

Differently, the study of the diffractive Drell-Yan reaction performed in [5j within the 
light-cone dipole description revealed importance of soft interactions with the partons spec- 
tators, which contributes on the same footing as hard perturbative ones, and strongly violate 
QCD factorization. 

Absorptive corrections are normally associated with the soft interactions between target 
and projectile, and they play an important role in diffractive hadron-hadron scattering. One 
can derive a Regge behavior of the diffractive cross section of heavy photon production in 
terms of the usual light-cone variables, 



i. 

pi ' i>- 



X 7 i — — , X 7 2 — — , K 1 -^) 



so that x 7 ix 72 = (M 2 + k^)/s and x 7 i — x~ j2 = x^p, where M, kp and x y p are the invariant 
mass, transverse momentum and Feynman xp of the heavy photon (di-lepton). n the limit 
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FIG. 1: The cross section of the diffr active DY process summed over all excitation channels at 
fixed effective mass Mx (left panel) corresponding to the Mueller graph in Regge picture (right 
panel) . 



of small x 7 i — > and large z p = vXlvX ~~ >* 1 the diffractive DY cross section is given by 
the Mueller graph shown in Fig. [TJ In this case, the end-point behavior is dictated by the 
following general result 

where ajp(t) is the Pomeron trajectory corresponding to the t-channel exchange, and e is 
equal to 1 or 1/2 for the Pomeron IP or Reggeon M exchange corresponding to 7* emission 
from sea or valence quarks, respectively (see Fig. CQ. 

As an alternative to the factorization based QCD approach, the dipole description of the 
QCD diffraction, was presented in Refs. [llj (see also Ref. [l2j])- It is based on the fact that 
dipoles of different transverse size r±_ interact with different cross sections cr(r±), leading to 
the single inelastic diffractive scattering with a cross section, which in the forward limit is 
given by [Tl[, 



dp 



167T 

P±=0 



where p± is the transverse momentum of the recoil proton, cr(r±) is the universal dipole- 
proton cross section, and operation (...) means averaging over the dipole separation. 

The color dipole description of Drell-Yan inclusive process first introduced in Ref. [l3j 
(see also Ref. |14|). treats the production of a heavy di-lepton like photon bremsstrahlung, 
rather than qq annihilation. Such a difference is a consequence of Lorentz non-invariance of 
the space-time description of the interaction, which varies with the reference frame. Only 
observables must be Lorentz- invariant. 

The dipole approach applied to diffractive Drell-Yan reaction in Ref. j5[, led to the QCD 
factorisation breaking, which manifests itself in specific features like a significant damping 
of the cross section at high y/s compared to the inclusive DY case. This is rather unusual, 
since a diffractive cross section, which is proportional to the dipole cross section squared, 
could be expected to rise with energy steeper than the total inclusive cross section, like it 
occurs in the diffractive DIS process. At the same time, the ratio of the DDY to DY cross 
sections was found in Ref. [EJ to rise with the hard scale, M 2 . This is also in variance with 



diffraction in DIS, which is associated with the soft interactions 15, 16 



The absorptive corrections affect differently the diagonal and off-diagonal terms in the 



hadronic current [171], in opposite directions, leading to an unavoidable breakdown of the 
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QCD factorisation in processes with off-diagonal contributions only. Namely, the absorptive 
corrections enhance the diagonal terms at larger y/s, whereas they strongly suppress the off- 
diagonal ones. In the diffractive DY process a new state, the heavy lepton pair, is produced, 
hence, the whole process is of entirely off-diagonal nature, whereas in the diffractive DIS 
contains both diagonal and off-diagonal contributions [ij]. This is the first reason why the 
QCD factorisation is broken in the DDY reaction. 

The second reason of the QCD factorisation breaking is more specific and concerns the 
interplay of soft and hard interactions in the DDY amplitude. In particular, this leads to 
the leading twist nature of the DDY process, whereas DDIS is of the higher twist jEf. Large 
and small size projectile fluctuations contribute to the diffractive DY process at the same 
footing, which further deepens the dramatic breakdown of the QCD factorisation in DDY. 
We will shortly discuss this issue below when presenting the numerical results. 



Quasieikonal model KMR for the so-called "enhanced" probability S en h (see e.g. Refs. [18 



19]), frequently used to describe the QCD factorisation breaking in diffractive processes, 



is not well justified in higher orders, whereas the color dipole approach considered here, 
correctly includes all diffraction excitations to all orders Q 1 . 

In this work, we investigate further the unusual features of the Drell-Yan diffraction 
[5| in the framework of the color dipole approach. We show that the unitarization effects 
can be correctly taken into account through eikonalization of the universal "bare" elastic 
dipole-target amplitude. This generalized dipole approach pretends to take into account 
the soft absorptive effects on the same footing with the hard dipole-target scattering. Such 
effects are included into the phenomenological partial elastic dipole amplitude fitted to data. 
This allows to predict the diffractive DY cross section completely in terms of the single 
parameterization of the dipole cross section known independently from the soft hadron 
scattering data. 

The paper is organized as follows. Section II contains derivation of the diffractive Drell- 
Yan amplitude in the dipole approach. Section III is devoted to a discussion of the unitarity 
corrections through the eikonalization of the elastic "bare" dipole-target scattering ampli- 
tude. In Section IV a short overview of different parameterizations of the elastic dipole-target 
scattering amplitude for small and large dipoles is given. The formulae for the single diffrac- 
tive Drell-Yan cross section are explicitly derived in Section V. Discussion of numerical 
results for the differential distributions and basic features of the diffractive DY is presented 
in Section VI. Finally, Section VII with the summary and conclusions closes our paper. 



II. DIFFRACTIVE DRELL-YAN AMPLITUDE IN THE DIPOLE APPROACH 

In the forward limit pt = 0, the photon radiation from a quark in inelastic collisions 
vanishes which was explained at the intuitive level, as well as demonstrated by a direct 
calculation of Feynman graphs in Ref. j20[. The same is also true for forward diffraction 
proceeding via gluon pair exchange (in the color singlet state) with no momentum transfer 
between the projectile quark and the proton target [5|. Disappearance of both inelastic 
and diffractive forward photon radiation happens due to the fact that if the electric charge 
gets no "kick", i.e. is not accelerated, no photon is radiated, provided that the radiation 
time considerably exceeds the duration time of interaction. This is dictated by the renown 



We are thankful to J. Bartels for pointing at this issue. 
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Landau- Pomeranchuk principle |2l| : radiation depends on the strength of the accumulated 
kick, rather than on its structure, if the time scale of the kick is shorter than the radiation 
time. 




a qq(n ~ r 2 ) Vqqin -r- 2 + af) 



FIG. 2: Leading order contribution to the diffractive Drell-Yan in the dipole-target collision. 



The non-Abelian case, QCD, is different: a quark can radiate gluons diffractively in the 
forward direction. This happens due to possibility of interaction between the radiated gluon 
and the tar get. Such a process, in particular, is very important for diffractive heavy flavor 
production [221 ] . 

Notice that the disappearance of Abelian radiation is only true for the diffractive scat- 
tering of a single quark off the target, and it does not hold for diffractive hadron-hadron 
scattering. As was demonstrated in Ref. js}, due to the internal transverse motion of the 
valence quarks inside the proton, which corresponds to finite large transverse separations be- 
tween them, the forward photon radiation does not vanish. This means that even at a hard 
scale the Abelian radiation is sensitive to the hadron size due to a dramatic break down of 



QCD factorization 23]. It was firstly found in Refs. |24l. |25| that factorization for diffractive 
Drell-Yan reaction fails due to the presence of spectator partons in the Pomeron. In Ref. j5[ 
it was demonstrated that factorization in Drell-Yan diffraction is even more broken due to 
presence of spectator partons in the colliding hadrons. 

As usual, we work in the rest frame of the target proton which remains intact after the 
collision. The hard part of the Drell-Yan process is given by the inelastic amplitude of 
7* radiation by a projectile quark (valence or sea) due to its interaction with the target 
through a gluon exchange as shown in Fig. [2j It consists of two terms corresponding to 
interaction of two different Fock states with the target - a bare quark before the photon 
emission \q) (s-channel diagram), and a quark accompanied by a Weizacker- Williams photon 
| g7*) (li-channel diagram). 

Let us consider first heavy photon bremsstrahlung by quark scattered off the proton 
target. We imply that the longitudinal momentum of the projectile quark cannot be changed 
significantly at high energies. In the high energy limit the corresponding s and w-channel 



contributions can be written as follows 13, 26 



~ -iZ q e a(l - a) „ 2f 2 + ^ ^ ( k 



^ zZ g ea]P % q Yl aXVl ] ***f(k±), ^ = (1-«)M 2 + «S (2.1) 
^7 a 2 (l± + k±y + r] 2 
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where k± = pj± + q± — p,ij_ is the transverse momentum of exchanged gluon, and l± = 
Pf± ~ (1 — a)q±/a is the transverse momentum of the final quark in a frame where z-axis is 
parallel to the photon momentum, and the amplitude for scattering of a quark off a nucleon 
in the rest frame of the nucleon reads 

Aw* (k±) ~ 2 P % ai V q {k±), A aaf (k ± ) ~ 2p° f 8 aaf V q {k ± ) 

Finally, we can switch to impact parameter space performing the Fourier transformation 
over l± and k± and write down the total amplitude for the photon radiation in inelastic 
quark-proton scattering 

M,(b,r) = -2ip° -Jii ' /1 ~ ° *%(a, r) ■ [v a (6) - V,(b + af) 



where af corresponds to the transverse separation between initial and final quark induced 
by the hard photon radiation, and ty^* q (a, r) is the light-cone wave function of the q — > 7*g 
transition in the mixed representation defined as follows 

d 2 L -r - u a {Pfh^Uaipf + q) 



/a I 



a 2 l\ + r/ 2 



The explicit expressions of the LC wave functions products for radiation of longitudinally 
(A = 0) and transversely (A = ±1) polarized photons are [HI, 14, 27 



2 

A=±l f/fj 



= Z l^{ m > 4K ° (m) K ( W2 ) + [1 + (1 - a) 2 ] ^^^K! ( m ) K x ( W2 ) 
^,(a,pi)*^ g (a,^) = ^<(A = 0)tt£ ff (a,pi) ei/ (A = 0)*£,(a, p 2 ) 

(J/CTj 

= Zj^M 2 (1 - a) 2 K ( m ) K ( W2 ) . 
Now let us turn to elastic dipole scattering as depicted in Fig. [2j It corresponds to forward 



scattering at small momentum transfers in the t-channel. Generally speaking, y^i — > Aqcd 
corresponds to the physical forward scattering limit since transverse momentum of a proton 
in the final state cannot be resolved to a better accuracy than its inverse size. 

In the leading order the elastic scattering amplitude is given by one-loop diagram with 
two it-channel gluon exchanges. Due to on-shell intermediate spectators, corresponding four- 
dimensional loop integral can be reduced to two-dimensional one over the transverse mo- 
mentum of one of the gluons 

2ilmF el (6 ± ) = J f^A(k ± )A(5 ± -k ± ), 5 ± <^\k ± \ (2.2) 

where A represents (inelastic) amplitude for one t-channel gluon exchange, and the last 
strong inequality guarantees that the proton target survives the scattering, hence, the elastic 
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nature of the process. Then the convolution theorem of Fourier analysis leads to the optical 
theorem 

lmF el (5) = J d 2 be-^ rb lmf el (b), 2iTm/ rf (&) = \A(b)\ 2 (2.3) 

which will be used below for eikonalization of multiple elastic amplitudes. 

Repeating calculations in this case we arrive at the qq dipole scattering amplitudes for s 
and it-channel photon emission, respectively, 



^ i 2 - 1 v *$ + ^ i 2 

C X CfCi 



x 

N, 



a 

^EE(l y ^-K(H f p + ar) 1 2 - I F 9 (6 + r; + ar) | 2 



X etc 



fCi 



where the last terms subtract the contributions from diagrams corresponding to the situation 
when none of the gluons couple to the same quark line with the hard photon. Then, implied 
the fact that all fields disappear at infinite separations, i.e. V q {b) — >■ when |6| — > oo, we 
have due to antisymmetry of the integrand 

/^"[ivy^r - \V,(b + r v + HI 2 "] - 0, n > 1 , Ifj , (2.4, 

such that these terms do not contribute to the final result. Using the optical theorem for 
the elastic amplitude 



2ilm/ ei (6,f p ) = \V q {b)-V q {b + r-)\\ 

X CfCi 



we can finally write 



Mff (6, f p , f, a) = -2ip° v 7 ^ ^ " *%(a, f) [21m / d (6, t>) - 21m / ei (&, f p + of)] (2.5) 



i.e. the amplitude of the diffractive radiation is proportional to the difference between elastic 
amplitudes of the two Fock components, with and without the photon radiation. When a 
quark fluctuates into the upper Fock quark-photon state with the transverse separation r, 
the final quark gets a transverse shift Ar = ar. Then the quark dipoles with different sizes 
in the \2q) and |2g7*) components interact differently, and their difference corresponds to 
the diffractive Drell-Yan process amplitude (12.51) . 



III. BREAKDOWN AND RESTORATION OF UNITARITY 

The elastic hadron-hadron scattering, which is observed in an experiment, is a compli- 
cated process, which can be composed to many elementary ( "bare" ) elastic scatterings which 
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are, in fact, the shadows of many inelastic interactions to be resummed to all orders. Such a 
resummation of the elementary scatterings, leads to unitarization of the elastic amplitude. 

By definition, an eigenstate of interaction cannot be diffractively excited, so can experi- 
ence only multiple elastic interactions. Correspondingly, its interaction cross section can be 
eikonalized, and this is not an approximation (like Glauber model for hadronic interaction), 
but is the exact result. In high-energy QCD the set of eigenstates are identified with the 
Fock states which can be treated as color dipoles. In the Drell-Yan reaction the lowest Fock 
state is an effective \qq) dipole [l3j]. Higher Fock states, like \qqg), etc. also contribute 
and their amplitudes should also be eikonalized. The approximation used here is to neglect 
those corrections. This is is justified for not very small fraction x 7 i and scale M 2 , where 
valence/sea quarks are dominated and the gluon contribution is rather small. 

In terms of the Regge theory, the elementary elastic dipole-target scattering corresponds 
to an exchange of the "bare" Pomeron. Assuming that this bare Pomeron is a Regge pole 
with the intercept above one (to justify the observed rising energy dependence of the cross 
section) one breaks down unitarity (Froissart bound) in the high energy limit. Eikonaliza- 
tion of this "bare" amplitude restores unitarity. The mentioned above higher Fock states 
correspond in Regge description to so called enhanced Regge graphs. Summing up all such 
graphs one arrives to the "effective" Pomeron, called Froissaron (28[, which satisfies the 
unitarity restrictions in both s-channel and t-channels. An additional justification to the 
mentioned above approximation neglecting higher Fock states, or enhanced graphs, comes 
from the observed smallness of the triple-Pomeron coupling. 

Having all this in mind, we can now easily generalize the expression (12.51) for the case of 
n consequent "bare" elastic scatterings of spectators as shown in Fig. [3J The amplitude in 
this case is given by 

TV. TV. 

The sign factor (— l) n is related to the phase of the bare Pomeron amplitude, which is nearly 
7r/2 due to the absorptive origin of the Pomeron, which is generated in the elastic amplitude 
by inelastic collisions through the unitarity relation. This can be also seen as a consequence 
of QCD as a non-Abelian theory [4]. Indeed, if it were an Abelian theory, the Born graph 
for the elastic amplitude would be one gluon exchange, i.e. the amplitude would be real. 
Since, however, QCD is non-abelian, the minimal number of exchanged gluons is two, i.e. 
the amplitude is imaginary. 

In (13.11) we followed neglected the contributions with the photon radiated by the quark be- 
tween subsequent elastic scatterings, in accordance with the Landau-Pomeranchuk principle 
mentioned above, which is at work if the radiation with a long coherence length, 

lc = 2^r y> L ' (3 - 2) 

where x 7 2 is defined in (II. 2p and L is the longitudinal distance covered by the interaction. 
In this case the radiation between the multiple interactions is suppressed by interferences 
(compare with gluon radiation in j29[), and the radiation spectrum depends only on the total 
momentum transfer, rather than on its multiple structure. Therefore multiple interactions 
do not lift the ban on a forward diffractive Abelian radiation, in spite of nonzero momenta 
transferred in each of the multiple scatterings. 
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Resummation of amplitudes (13. ip to all rescattering orders results in the total amplitude, 
which has the following general form 



MUb,r p ,r,a) = 2ip»V4K 



\/l — a 



e -2Im/ e i {b,f p ) _ e -2Im/ ei {b,r v +a>r) 



(3.3) 



Note, in this derivation we implied that relation (12. 4p holds in each particular order n, thus 
the subtraction terms do not contribute to the final expression (13. 3p . 




FIG. 3: Total eikonalized amplitude of the diffractive Drell-Yan process, which effectively includes 
the gap survival effects. 



Let us now write the total hadronic amplitude of the diffractive Drell-Yan process as |5| 

Aif = 4? + A® + A® , (3.4) 
where each term corresponds to 7* radiation by one of the valence quarks, in particular, 

x ^ i (f u f 2 ,f 3 ;x qi ,x q2 ,x q3 )^f* f ('r 1 + af,f 2 ,f 3 ;x qi - x jX , x q2 ,x q3 ) 

x Mq q (b, n - f 2 , r, a) + Mq q (&, n - f 3 , f, «)] e ^' a V^- g (3.5) 

Here, A 7 = L,T; l± = 5± — q±/at (z-axis is directed along initial proton momentum); the 
hard photon with virtuality q 2 = M 2 3> m 2 , transverse q± and fractional longitudinal x 7 x 
momenta is emitted from the first valence quark with impact parameter f\ (see Fig. [3]), 
other two valence quarks in the proton have impact parameters r 2 and r 3 , respectively; 
r is transverse separation between the photon and the radiating quark; a = x 1 \/x qi is 
the fraction of longitudinal momenta taken away by the photon from the radiating quark; 
M qq are the Fourier-transformed amplitudes for the elastic quark dipole scattering off the 
proton target accompanied by the hard L, T-polarized photon emission calculated above in 
Eq. (I3.3P ; are the light-cone wave functions of the 3q systems in the initial and final 
state, respectively. In Eq. (I3.5P we implicitly assumed that exchanges t-channel gluons all 
together take a negligibly small longitudinal momentum compared to the collisions energy 
a/s and, hence, corrections to quark momenta due to gluon couplings are neglected in the 
wave functions. 

Eikonalization of the universal elastic dipole-target scattering amplitude performed in 
Eq. (I3.3P incorporates all soft and hard interactions between qq dipole and the target on 
the same footing as the parameterization for the amplitude on a free proton, f e i(b, r), fitted 
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to data on DIS and soft hadron scattering. This provides an alternative to the conven- 
tional treatment of the gap survival effects in terms of the single suppression factor in the 
cross section cjddy = K ■ Uddy^ the so-called gap survival factor, estimated in the eikonal 
approximation 



30 



K _, 1 oZ(s) , 1 WZ(s)? 

vr BST(a) + 2B™{s) + (An) 2 E%(s)[B%(s) + B p J(s)} ^ 

where the energy-dependent elastic slope is B p e f(s) = B® L + 2a'jp ln(s/so) with = 
7.5 GeV" 2 , s = 1 GeV 2 . The slope of single-diffractive DY cross section can be estimated 
as, B^ ~ ( r ch)/3 + 2a'jp\n(s / so) , where the proton mean charge radius squared (r 2 h ) = 0.8 
fm 2 . More elaborated models for the gap survival factor incorporating a part of the Gribov 



corrections (see e.g. Refs. |31|, |32|) predict similar suppression factors, and one can easily 
replace the K-f actor (13. 6 j) by a preferable one. 

In order to demonstrate that the eikonalization procedure (13. 3p correctly takes into ac- 
count soft gap survival effects, we checked that the ratio between diffractive Drell-Yan cross 
sections with eikonalized (I3.3P and non-eikonalized (12. 5 j) diffractive amplitude leads to a 
suppression factor which is very close numerically to the standard gap survival factor K 
defined independently from Eq. (13. 6p . 



a l 

^noneik 
DDY 



K (3.7) 



and such a relation between them holds at different energies. And this ratio does not depend 
on Xyi and M 2 . This confirms our statement we have made in the beginning of this Section: 
the eikonalization of the elastic dipole-target amplitude correctly incorporates the unitarity 
corrections. 

As was already been said above, the parameterizations for the real elastic dipole-target 
scattering fitted to experimental data must already effectively contain the absorptive and 
all-order QCD corrections, so it should not be eikonalized when used in explicit calculations, 
and no iT-factors are to be applied to the cross section in this case. 

Let us now shortly discuss different parameterizations for the elastic dipole-target scatter- 
ing amplitude f e i(b,r) corresponding to the scattering of small (r r p ) and large (r ~ r p ) 
dipoles off the proton target, known from the data fits. 



IV. ELASTIC DIPOLE-TARGET SCATTERING AMPLITUDE 



It is well-known that small x (large Q 2 ) regime corresponding to a scattering of small 
dipoles with r — > is well described by the popular Golec-Biernat-Wuestoff (GBW) 
parametrization of the dipole cross section (33|. The elastic 6-dependent amplitude in this 
case is x-dependent and has a form 34] (see also Ref. (3o| ) 



lmf^ w (b,r,x,x q ) 



2 exp 



8ttB(x) 



exp 



[b + f(l -Xg)f 

2B(x) 



+ exp 



[b + fxg] 2 
2B(x) 



mix) 



[b + r(l/2-x q )]' 
2B(x) 



B(x) = R%(x) + R 2 (x)/8 



(4.1) 
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where parameters are fitted to DIS data at small x [33[ <To = 23.03 mb, Ro{x) = 
0.4 fm x (x/xo) ' 144 with xq = 3.04 x 10~ 4 , x ~ Q 2 /s is the Bjorken variable, x q is the 
quark longitudinal quark fraction in the dipole defined in Eq. (15. 9p . R 2 N (x) in the limit 
r — can be defined through the slope of elastic electroproduction of p-mesons measured 
at HERA as 

R%(x) = B rp ^ pp {x, Q 2 > 1 GeV 2 ) - R 2 {x)/4, B rp ^ pp {x, Q 2 > lGeV 2 ) ~ 5 GeV" 2 
. Amplitude (14. ip correctly reproduces the dipole cross section [33| 



2 J d 2 b hnf e i(r, b, x, x q ) = a qq (r, x) 



but contains more information, because it is sensitive to the color dipole orientation within 
the phenomenological saturation model, which includes contributions from higher order 
perturbative corrections as well as non-perturbative effects contained in DIS data. 

Notice that the simple GBW parameterization (14. ip has some restrictions. In particu- 
lar, the non-integrated gluon distribution exhibits no power-law tails in momentum space 
in contradiction with QCD. Moreover, it does not match the DGLAP evolution at large 
values of Q 2 . Therefore, one should be cautious applying this model at very high transverse 
momenta accessible at the energies of LHC j35| . 

For soft scattering (moderate and small Q 2 ) corresponding to large dipoles 
the cm. energy squared s, rather than Bjorken Xj_js the proper variable. So, for the soft 
processes one can switch from x- to s-dependence [20|, keeping the same functional form of 
the dipole amplitude (14. ip and adjusting the parameters to observables in soft reactions as 



35-37 



Imf% BW (b,r,x,x q ) lmf* ST (b,f,s,x q ), (4.2) 

R (x) ->■ Rq(s) = 0.88fm(s /s) - 14 , 

R 2 N (x) R 2 N (s) = BZ(a) - \r 2 (s) - \{r 2 ch )^ 



(T -> (r (s) =aZ(s)(l + 



8(r c \). 



where the pion-proton total cross section is parameterized as [38j cr^(s) = 23.6(s/so) ' 08 
mb, s = 1000 GeV 2 , the mean pion radius squared is j39| (r 2 ^^ = 0.44 fm 2 , and the 
Regge parametrization of the elastic slope B^(s) = B + 2a'jp\n(s/ /j 2 ), with B = 6 GeV -2 , 
a'jp = 0.25 GeV -2 , and /i 2 = 1 GeV 2 can be used. We shall refer to Im/? ST (6, r, s, x q ) 
below as the Kopeliovich-Schafer-Tarasov (KST) parametrization of the elastic dipole-proton 
amplitude (37j . 

Two models for the elastic quark dipole-target scattering amplitude f e i(b,f), given above 
by Eqs. (14.11) and (14.21) . are valid for small r ^ r p and large r ~ r p dipole (r p is the mean 
proton size), respectively. Since the diffractive DY cross section is primarily sensitive to the 
large transverse separations ~ r p between target and projectile implied by the forward limit 
t — > 0, then the KST parameterization for the dipole cross section should be used, at least, 
in the leading order calculation. For completeness, we will compare the DDY cross sections 
calculated with GBW and KST parameterizations below when discussing the numerical 
results. 
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V. SINGLE DIFFRACTIVE DRELL-YAN CROSS SECTION 



The differential cross section for the single diffractive di-lepton production in the target 
rest frame reads 



3 



f n=l -/ \—/ ~\Si,n; 



in a(l-a)(27T) 5 8(p° ^ 2 



x ^AKW^^, A 7 = L,T (5.1) 

where prefactors provide averaging over colors and helicities of exchanged t-channel gluons, 
p° n is the energy of the radiating rath quark in the initial state, n — 1, 3; a em = e 2 / (4m) = 
1/137 is the electromagnetic coupling constant. The second line in Eq. fl5.ip describes decay 
of 7* into the leptonic pair l + l~ into solid angle dQ = d<fid cos 9, and the standard leptonic 
tensor is given by 

IT = 4(pf +P r- + tfi+vl - gTiPi+ ■ Pi-)) ■ 

At the moment, we are not interested in lepton polarizations and their angular distributions, 
so for the sake of simplicity we integrate out the cross section (15.11) over the solid angle of 
the lepton pair. We keep in the cross section only diagonal in the photon polarization 
A 7 = L,T terms (non-diagonal ones drop out after integration over leptonic azimuthal 
angle <fi). Integrating the diffractive differential DY cross section over the photon transverse 
momentum q± we get 



(5.2) 



d 4 (?L,T{pp pl + l X) _ a em d 3 a L , T (pp -»• P1*X) 
dlnadM 2 d 2 5 ± ~ 3vrM 2 dlnad 2 5 ± 

Then applying the completeness relation 

5^#/(ri + ar,r 2 ,r 3 ;x qi ,x q2 ,x q3 )V* f ( r x + ocr , r 2 1 t 3 ; x qi , x q2 , x q3 ) 
f 

3 

= S(n - r\ + a(r- f'))8(f 2 - f 2 )5(f 3 - f 3 ) JJ<K%- - x' q .) (5.3) 

i=i 

we get the diffractive 7* production cross section in the following differential form 

A ^ ^ "' ^ d 2 rid 2 r 2 d 2 r 3 d 2 r d 2 bd 2 V dx qi dx„ 2 dx q 



dlnad 2 S ± 64vr 2 J L 1 6 ^ qi ^ q2 ^ q3 

ot,r)\ |^j(ri, r 2 , r 3 ; x qi , x q2 , x qz )\ 
x A(fi, f 2 , f 3 ; b; r, a)A(f 1} f 2 , f 3 ; b'; r, a) e lS ^< b - b ') (5.4) 



where \E , 7 . g = \l/ 7 » 9 /Z 9 , and 



A = -21m /| ST (6, ft - f 2 ) + 21m f% ST (b, f x -f 2 + af) 

-21m / e f T (6, n - f 3 ) + 21m / e f T (6, f x - f 3 + af) , (5.5) 
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where the KST parameterization (14. 2 j) fitted to the soft data and hence valid at \fi — fj\ 
b, i ^ j is used. Finally, going over to the forward limit 5± = we get 



d 3 a Xy (pp ^ m*X) 



din a dS , 



s± =o 



64vr 



d 2 rid 2 r 2 d 2 r 3 d 2 r dx qi dx q2 dx q3 



7 3 



gu M qii -^qz, 



-i 2 



rf 2 6A(fi,f 2 ,f 3 ;6;f,a) 



(5.6) 



We see that normalization of the cross section agrees with the original result of Ref . [5] • The 
total diffractive cross section is then given by 



da(pp — > P7*X) 
din a 



d 3 a(pp — > P7*X) 



din a dd 2 ^ 



(5.7) 



where B^(s) is the diffractive slope similar to the one measured in diffractive DIS. 

The next step is to introduce the proton wave function assuming the Gaussian shape for 
the quark distributions in the proton as 



I 5^2;^; X q 1 j X q 2 1 X q3 , 



2 + a/6 
ir 2 ab 



exp 



r 2 + r 2 



33 1 



x 5(fi + r 2 + r 3 )S(l 



x. 



qi 



33 / 



(5.8) 



where a = (r? g ) and 6 = (R 2 ) are the diquark mean radius squared and the quark mean 



distance from the diquark squared, respectively. In this work, we will use the simplest case 
of symmetric valence quarks distribution assuming that r qq = 
Then valence quark distribution in the proton is given by 



R q = 0.85 fm. 



/ 



dx q2 dx q3 p(x qi , x q 2 , x qs ) p q ^(x 



31 - 



where we integrated out the longitudinal fractions of the diquark in the proton. Generaliza- 
tion of the three-body proton wave function (15.81) including different quark and antiquark 
flavors leads to the proton structure function [40| 



1 



Z q[Pq( X q) + Pq( X q)} = — F i( X q) 



,1, 



a 



(5.9) 



In the numerical analysis below, in order to estimate the theoretical uncertainty in the 
diffractive DY process we will use a few different parameterizations for the proton structure 
function F 2 , widely used in the literature. 



VI. NUMERICAL RESULTS 

Let us now turn to discussion of the numerical results. We start from the comparison 
of the differential cross sections for single diffractive and inclusive Drell-Yan processes. In 
Fig. H]the ratio of the diffractive to inclusive DY cross sections is plotted as a function of di- 
lepton invariant mass squared M 2 (left panel) and photon fractional light-cone momentum 
x 7 i (right panel) at different energies. In the left panel, the curves are given for fixed 
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Xji = 0.5 (solid lines) and x 7 i = 0.9 (dashed lines). In the right panel, the curves are given 
for fixed M 2 = 50 GeV 2 (solid lines) and M 2 = 500 GeV 2 for y/H = 14 TeV and 500 GeV, and 
M 2 = 200 GeV 2 at y/s = 40 GeV (dashed lines). The pairs of solid/dashed curves in the both 
panels correspond to \/s = 40 GeV, 500 GeV and 14 TeV from top to bottom, respectively. 
Here we used the KST parameterization for the dipole-target scattering amplitude (35-33 
and F2 parameterization by Cudell and Soyez jUJ are used here and below unless otherwise is 
specified. Unitarity corrections are included by default into the employed phenomenological 
partial elastic dipole amplitude fitted to data (see above). In this calculation we consider the 
unpolarized case summing up the contributions of longitudinal and transverse parts both in 
the diffractive and inclusive cross sections. 



do"DDY / dx r i dM 2 
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FIG. 4: The ratio of the diffractive to inclusive Drell-Yan cross sections as function of the lepton- 
pair invariant mass squared M 2 (left panel) and photon fraction x~\ (right panel) at different 
energies. In the left panel, the curves are given for fixed x 7 i = 0.5 (solid lines) and x-yi = 0.9 
(dashed lines). In the right panel, the curves are given for fixed M 2 = 50 GeV 2 (solid lines) and 
M 2 = 500 GeV 2 for = 14 TeV and 500 GeV, and M 2 = 200 GeV 2 at = 40 GeV (dashed 
lines). The pairs of solid/dashed curves in the both panels correspond to yfs = 40 GeV, 500 GeV 
and 14 TeV from top to bottom, respectively. 



As seen from Fig. HJ the DDY-to-DY cross section ratio is falling with energy. However, 
naively one could expect basing on QCD factorisation, that the DDY cross section, which 
is proportional to the dipole cross section squared, should rise with energy steeper than the 
total inclusive cross section. At the same time, the ratio rises with the hard scale of the 
process, M 2 . This also looks counterintuitive, since diffraction is usually associated with soft 



interactions [13j . These effects are different from ones emerging in Regge factorisation-based 
calculations, where we observe a slow rise of the DDY-to-DY cross section ratio with c.m.s. 
energy and its practical independence on the hard scale of the process M 2 

In order to understand such an interesting shape of the energy and hard scale dependence 
of the DDY-to-DY cross section ratio obtained in the color dipole approach, let us look at 
the amplitude of the DDY process, which is proportional to the difference between the dipole 



cross sections of the Fock states with and without the hard photon emission [5J, i.e. 

2a<7o 



M DDY ~ *{R) - a{R - of) = ^-f R2/R2 ° ix) (r-R)+ h.o. (6.1) 



assuming the simplest GBW slope for the dipole cross section, and the hardness of the 
emitted photon implies r ~ 1/M Rq(x). We see now that the diffractive DY amplitude 
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is linear in r, so the diffractive cross section turns out to be a quadratic function of r, which 
is different from e.g. the diffractive DIS process where the cross section is proportional to 
r 4 and is dominated by soft fluctuations (see e.g. Refs. [3, EEf). Since the diffractive DY 
cross section is proportional to r 2 , then soft and hard interactions contribute on the same 
footing [5|, which is one of the basic sources of the QCD factorisation breaking in diffractive 
DY process. 

As was demonstrated in Ref . j5[ , all the energy and scale dependence of the DDY-to-DY 
cross section ratio comes via the x-dependent factor, 

^ oc 1 e -a*/*8<») (6.2) 
a DY R$(x) 

where x = M 2 /x 7l s. In our case, R%(x) < 2R 2 , so the factor in Eq. (16.21) rises with Ro(x), 
i.e. with x. This is the main reason why the ratio shown in Fig. H] decreases with energy, 
but increases with the hard scale M 2 . Also, the falling energy behavior is partly due to the 
rise with energy of the absorptive corrections 



do^DDv/dXyidM 2 

dcr T DDY / dx y i dM 2 
0.200 
0.150 

0.100 ^ 
0.070 
0.050 

0.030 

0.020 
0.015 

0.010 



M 2 , GeV 2 



do" L DDY/dx yl dM 2 

do" F DDY / dx yl dM 2 
0.200 
0.150 

0.100 
0.070 
0.050 ■ 



M 2 =50 GeV 2 





0.30 



0.50 



0.70 



1.00 



FIG. 5: The ratio of the longitudinal (L) to transverse (T) photon polarization contributions to 
the diffractive Drell-Yan cross section as function of the lepton-pair invariant mass squared M 2 
(left panel) and photon fraction x 7 i (right panel) at different energies: y/s = 40 GeV (dash-dotted 
lines), y/s = 500 GeV (solid lines) and y/s = 14 TeV (dashed lines). In the left panel, the curves 
are given for fixed x 7 i = 0.5 (upper three curves) and x 7 i = 0.9 (lower three curves). In the right 
panel, the curves are given for fixed M 2 = 50 GeV 2 . 



In Fig. [5] we show the relative contribution of the longitudinal (L) to transverse (T) 
photon polarization to the diffractive Drell-Yan cross section. The ratio ol/ot is presented 
as function of lepton-pair invariant mass squared M 2 (left panel) and photon fraction x 7 i 
(right panel) at different energies: y/s = 40 GeV (dash-dotted lines), y/s = 500 GeV (solid 
lines) and y/s = 14 TeV (dashed lines). In the left panel, the curves are given for fixed 
x 7 i = 0.5 (upper three curves) and x 7 i = 0.9 (lower three curves). In the right panel, 
the curves are given for fixed M 2 = 50 GeV 2 . We see that the diffractive DY process is 
always dominated by radiation of transversely polarized lepton pairs. The ratio ul/ctt is 
only slightly dependent on M 2 , and there is no any significant energy dependence. The 
longitudinal photon polarization amounts to about 10 % at x 7 i = 0.5 and then steeply falls 
down at large x 7 i — > 1. Such a behavior turns out to be similar to that for inclusive DY 
process 
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V7=500GeV 2 V7=500GeV 

d<r DDY / dx 7 ,dM 2 , fb dcr DDY/ dx rl dM 2 , fb 

1000 r 




FIG. 6: Diffractive Drell-Yan cross section (in fb) as function of the lepton-pair invariant mass 
squared M 2 (left panel) and photon fraction x*yi (right panel) at the RHIC II c.m.s. energy 
y/s = 500 GeV. In the left panel, the curves are given for fixed x 7 i = 0.5 with GBW (long- 
dashed line) and KST (solid line) parameterizations, and for x~i = 0.9 with GBW (dash-dotted 
line) and KST (dashed line) parameterizations. In the right panel, the curves are given for fixed 
M 2 = 50 GeV 2 with GBW (long-dashed line) and KST (solid line) parameterizations, and for 
M 2 = 500 GeV 2 with GBW (dash-dotted line) and KST (dashed line) parameterizations. 



We also compare predictions for the diffractive DY cross section for different parameter- 
izations for elastic dipole-target scattering amplitude corresponding to scattering of small 
(GBW given by Eq. (14. ip ) and large (KST given by Eq. (14.21) ) dipoles. As an example, 
in Fig. [6] we present the diffractive Drell-Yan cross section as function of the lepton-pair 
invariant mass squared M 2 (left panel) and photon fraction x 7 i (right panel) at the RHIC II 
c.m.s. energy yfs = 500 GeV. We notice that the GBW parameterization leads to roughly 
a factor of two smaller cross section than the one obtained with the KST parameteriza- 
tion, however, both of them exhibit basically the same x 7 i and M 2 shapes. It means that 
the evolution of the dipole size can only affect the overall normalization of the DDY cross 
section. Since arguments in the elastic amplitude f e i in Eq. (I3.3p . the impact distance be- 
tween the target and the projectile b and the transverse distance between projectile quarks 
r p ~ l^i — fj\, i j, are of the same order and given at the soft hadronic scale, then the 
use of KST parameterization fitted to the soft hadron scattering data data is justified in the 
case of diffractive DY. 

In order to illustrate the intrinsic theoretical uncertainties in our DDY cross section cal- 
culations, in Fig. [7] we show the diffractive DY cross section as function of the lepton-pair 
invariant mass squared M 2 (left panel) and photon fraction x 7 i (right panel) for differ- 
ent parameterizations of the proton structure function F 2 entering the DDY cross section 
through Eq. (15. 9p . In this figure we represent results with four distinct cases widely used 
in the literature: Regge parameterization by Cudell and Soyez [ijj (solid line), old SMC 
parameterization |42| with infrared freezing at two different scales Qa, = 0.1 (dashed line) 
and 0.5 GeV 2 (long-dashed line) and recent GJR parameterization |43| (dash-dotted line). 
We see that the diffractive DY cross section is sensitive to various F 2 parameterizations, 
especially at relatively large x 7 i — >■ 1, which is reflected in quite noticeable uncertainties in 
our calculation. This opens up a promising opportunity to use the diffractive DY reaction 
as a direct probe of the proton structure function at rather large j a. 

The Regge parameterization 41] is presumably constructed in the soft region where the 
most of the contribution to the DDY process comes from, and it leads to a rather regular 
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FIG. 7: Diffractive Drell-Yan cross section (in fb) as function of the lepton-pair invariant mass 
squared M 2 (left panel) with fixed = 0.5 and photon fraction x-y\ with fixed M 2 = 50 GeV 2 
(right panel) for different par ameterizations of the proton structure function Fi : Regge parame- 
terization by Cudell et al 41] (solid line), SMC parameterization [i^] with infrared freezing at two 
different scales Q\ = 0.1 (dashed line) and 0.5 GeV 2 (long-dashed line) and GJR parameterization 



431 ] (dash-dotted line). 



and stable behavior of the cross section at x 7 i — > 1. Old SMC parameterization is 
obviously inapplicable to the DDY calculations, because it leads to significant uncertainties 
with respect to the lower freezing scale Qo variations indicating at a strong sensitivity to 
the non-perturbative low-Q region where the proton structure function F2 is unknown to a 
large extent. 



VII. CONCLUSION AND OUTLOOK 

In this work, we have investigated in detail the QCD factorisation breaking effects in 
the diffractive Drell-Yan process within the framework of the color dipole approach. Such 
effects lead to quite different properties of the corresponding observables with respect to 
QCD factorisation-based calculations. 

A quark cannot diffractively radiate a photon in the forward direction, whereas a hadron 
can due to the presence of transverse motion of spectator quarks in the projectile hadron. 
For this reason, the diffractive DY cross section depends on the hadronic size breaking the 
QCD factorisation. 

This leads to the physical picture where hard and soft interactions are equally important 
for DY diffraction, and their relative contributions are independent of the hard scale, like 
in the inclusive DY process. This is a result of the specific property of DY diffraction: 
its cross section is a linear, rather than quadratic function of the dipole cross section. On 
the contrary, diffractive DIS is predominantly a soft process, because its cross section is 
proportional to the dipole cross section squared. 

Contrary to what follows from the calculations based on QCD factorisation, the ratio 
of the diffractive to inclusive cross sections falls with energy, but rises with the di-lepton 
effective mass M. This happens due to the saturated behavior of the dipole cross section 
which levels off at large separations. All these properties are different from those in the 
diffractive DIS process, where QCD factorisation is exact. 

In addition, we made predictions for the differential (in photon fractional momentum x 7l 
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and di-lepton invariant mass squared M 2 ) cross sections for the diffractive DY process at the 
energies of RHIC (500 GeV) and LHC (14 TeV). The transverse photon polarisation gives 
the dominant contribution to the DDY cross section. The ratio <Jl/&t is almost independent 
on c.m.s. energy y/s and only weakly depends on the hard scale M 2 . 

Finally, we propose an alternative treatment of the absorptive effects describing the sub- 
sequent soft interactions of projectile quarks off the proton target by multiple elastic dipole- 
target scatterings. Such an idea leads to eikonalization of the "bare" elastic dipole amplitude 
in the DDY amplitude and, ultimately, to a description of the gap survival effects on the 
same footing with the leading-order diffractive DY subprocess in the framework of the color 
dipole approach without introducing the gap survival factor K. These corrections are in- 
cluded by default into the employed phenomenological partial elastic dipole amplitude fitted 
to data. 

The main features of the diffractive Drell-Yan reaction described above are valid for other 
diffractive Abelian processes, like production of direct photons, Higgsstrahlung, radiation of 
Z and W bosons. A detailed analysis of these processes in the framework of the color dipole 
approach is planned for a forthcoming study. 
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